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1 The experimental progress 

The experimental progress made in the 1990's is remarkable in many ways. With 
the right experiments we can reveal the fine details of high energy phenomena. On 
the eve of LEP shutdown we known that physics at the microscopic scale has a basis 
as rationale as chemistry, everything moves, not by magic, but because there is a 
mechanism at work. Everything is encouraging for the major unsolved aspects of the 
breaking of electroweak gauge symmetry. 

The blue band, giving the theoretical uncertainty in the ^-distribution as a func- 
tion of the Higgs boson mass is, perhaps, the most celebrated figure of the whole LEP 
era. It is interesting to examine the environment in which the technology for produc- 
ing this figure has been developed. One way is through the Kiihn tower. According 
to this dantescan representation we find at the top of the tower a group of highly- 
specialized experts: QED, QCD etc. Next to them we will find the Smugglers: those 
who assemble a code mooncursing it. At the lower level the Users are represented: 
the experimental community at large. 

The Rosetta stone that keeps this world running is based on some sort of universal 
language that has been created to allow for a meaningful communication, the Pseudo- 
Observable language. 

1.1 The language of Pseudo-Observables 

The PO are related to measured cross- sections and asymmetries by some deconvo- 
lution or unfolding procedure and the concept itself of pseudo-observability is rather 
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difficult to define. One way to introduce it is to say that the experiments measure 
some primordial (basically cross-sections and thereby asymmetries also) quantities 
which are then reduced to secondary quantities under some set of specific assump- 
tions. Within these assumptions the secondary quantities, the PO, also deserve the 
label of observability. 

A quick example of this language is as follows. The explicit formulae for the Zff 
vertex are always written starting from a Born-like form of a pre-factor x fermionic 
current, where the Born parameters are promoted to effective, scale-dependent pa- 
rameters, 

fa» [{if +ia L ) 1+ -2Q f Kf z s 2 + ia Q ] = 7 „ {Q{ + Q{ 7s ), (1) 

where 7+ = 1 + 75 and a Q L are the SM imaginary parts. 

By definition, the total and partial widths of the Z boson include also QED and 
QCD corrections. The partial decay width is therefore described by the following 
expression: 

r/ = r(z-/7)=4 C /. r (\gf v fRf v + \gf A fR{) + A EW/QCU , (2) 

where cf. = 1 or 3 for leptons or quarks (/ = I, q), and i?£ and W A describe the final 
state QED and QCD corrections and take into account the fermion mass mj. The 
last term, 

a =r (2) _^ir (1) ($\ 

"EW/QCD *■ EW/QCD _ EW ' \ / 
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accounts for the non-factorizable corrections. The standard partial width, r , is 

G ilf 3 

' ' " 82.945(7) MeV. (4) 



The peak hadronic and leptonic cross-sections are defined by 

z z z z 

where V z is the total decay width of the Z boson, i.e, the sum of all partial decay 
widths. The effective electroweak mixing angles ( effective sinuses) are always defined 
by 

Re a? 

4|Q / |sin 2 0i= i__t = i_^ > gf=BBgf, g{ = Re Q{ . (6) 
The forward-backward asymmetry A pg is defined via 

^fb = Z L ^T > a T = a F +a B , (7) 

F + (T B 

where a F and a B are the cross-sections for forward and backward scattering, respec- 
tively. 



2 A new era 

Admittedly, we are entering a new era: we have been moving from a mediaeval 
architecture of radiative corrections to the renaissance of LEP period and now one 
should not indulge in deciphering the intimate meaning of the last digit in the last 
fit of the LEPEWWG. 

A new frontier is at the horizon: most likely it is goodbye to the one man show. 
Running a new RC project will be a little like running an experiment. 

At the eve of LEP shutdown it is of some importance to summarize the present 
status of high precision physics |J. For e + e~ — > ff all one-loop terms are known, 
including re-summation of leading terms. At the two-loop level leading and next- 
to-leading terms have been computed and included in codes like TOPAZO @^J] and 
ZFITTER ||. For realistic observables initial state QED radiation is included via 
the structure function method, or equivalent ones. Final state QED is also available 
as well as the interference between initial and final states 0. Fine points in QED 
for 2^2 are as follows. For s-channel all the O (a 2 L n ) , n = 0,1,2 terms are 
known from explicit calculations, the leading O (a 3 L 3 ) is also available and they are 
important for the studies of the Z lineshape. 
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Differences and uncertainties amount to at most ±0.1 MeV on M z and T z and 
±0.01% on al (MIZA, TOPAZ0 and ZFITTER) For non-annihilation processes 
(Bhabha) both structure-function and parton-shower methods have been analyzed 
and the uncertainty is estimated to be 0.061% from BHLUMI [|7|]. Certainly, full 
two-loop electroweak corrections are needed for GigaZ (10 9 Z events) with a quest for 
a fast numerical evaluation of the relevant diagrams. 

For e + e~ — > 4fermions all tree-level processes are available and O (a) electroweak 
corrections are known only for the W^M^-signal and in double-pole approximation 
(DPA) § and |J. e + e~ — > 4f + 7 in Born approximation is also available for all 
processes [[l(|. Dominant radiative corrections are inserted in single- W production 



through the so-called Fermion-Loop scheme [ |TT| - |l~3| 



Fine points in QED for 2 — > 4 are as follows For e + e~ — > WW — > 4f DPA 
gives the answer but, for a generic process e + e~ — > 4f QED radiation is included by 
using s-channel structure-functions, i.e. in leading-log approximation. The latter 
are strictly applicable only if ISR can be separated unambiguously. Otherwise their 
implementation may lead to an excess of radiation. Preliminar investigations towards 
non-s SF by GRACE and by SWAP [ID| gives an indication on how to implement 
the bulk of the non- annihilation effect but still represent ad hoc solutions. These 
methods, which are essentially based on a matching with the soft photon emission, 
still contain an ambiguity on the energy scale selection with consequences on the 
predicted observables. 



3 Realistic observables 

How well do we know realistic observables in e + e~ — > 2f? Pretty well around 
the Z peak, reasonably well up to LEP 2 energies, but fermion-pair corrections (PP) 
are not really under control. Let us describe the structure of PP by starting with a 
simple case, e + e~ PP-corrections to e + e~ — > bb. PP is more than photon conversion, 
we will have the following sets of diagrams: Multi-Peripheral, Initial State Singlet 
(ISS), Initial State Non-Singlet (ISNS) and Final State (FS). 

Note that we include both 7 and Z exchange. On top of real pair production 
one has to include virtual e + e~ pairs. Let us compare as a typical example the 
sum of virtual and real pair corrections for primary hadrons and primary muons at 
y^s) = 189 GeV. 

The comparison between GENTLE, ZFITTER, and TOPAZ0 for the diagram- 
based definition reveals maximum differences of 1.7(1.5) per mill for inclusive hadrons 
(muons) and 0.2(0.4) per mill for high s' hadrons (muons). 

However, PP are to be completed expressly at higher energies and above vector 
boson thresholds. The first evolutionary development is: a 4f MC has to be fully 
interfaced with a 2f MC. 
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4 QED radiation 



There are several 4f processes, those with t-channel photons that are not domi- 
nated by annihilation. Examples are: single-VT production and two-photon processes. 
How to include the bulk of QED radiative corrections? The second evolutionary de- 
velopment: is multi-photon radiation a one-scale or a multi-scale convolution phe- 
nomenon? 

a ( p+p _ _> { qi } + QED) = f [] dx x D(x x , ?) a (x +P+ x^ -> {q,}) (8) 

In the above equation the question mark means that the corresponding scale has to be 
guessed. We need to understand how the standard SF-method is related to the exact 
YFS exponentiation. In the standard YFS treatment of multiple photon emission we 
have 

<rh> + +p_- E <H + E k j) ~ fdPS q \M \ 2 e( P+ + p _-]T % V (9) 

\ 1=1,2/ j=l,n J V i / 

where E is the spectral function defined through the usual eikonal factor: 
E(K) = j^- 4 J d 4 xexp(tK-x)E(x), 

E(x) = expj^ J d 4 ke lk - x 5 + (k 2 ) \ f{k) | 2 | (10) 

At this point we choose an alternative procedure were we do not separate the soft 
component from the hard one and compute some exact result valid for an arbitrary 
number of dimensions n and for on-shell photons, i.e. k 2 = 0, 



-/ 



Pi ■ kpj ■ k 



2. 



In dimensional- regularization one has the following result, valid \/x 

I(x) = -7r Pj [ 1 ^^ + 21n2-lnx 2 -ein^, (12) 

where we have defined a variable £ as the ratio £ = \xo\/r, with an infinitesimal 
imaginary part attributed to xq, i.e. xq — > xq + iS (5 — > 0+). Furthermore, P is 
the linear combination P = pj + (ppi — Pj) u, where we have defined p to satisfy 
(ppi — pj) 2 = 0, and xq, r are rewritten in covariant form as follows: 

'Ob Q O O , i 

x ° = ^^JTp^ ' r =x o + x - ( 13 ) 
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The last integral shows the infrared pole - and a collection of Li 2 -functions. There- 

e 

fore, E(K) is not available in close form. The scheme that we want to propose defines 



a coplanar approximation [14] to the exact spectral function, 



def 2 _ t 1 ^PijPi tc def 2 1 



3 3 pi-PiiPi Pi 3 



J^cp = In <j e" AlR Pi ' XPj ' X I A m = - + constants. (14) 
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Within the coplanar approximation we have 

i r e - AiR 



£ipair<ij>/^\ 
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X 



(2 7r) 2 [mjmjj ^(aAy) 

/ dcTrfcr' (ctct') J " <$ 4 (op» + O-'ft' - K) . (15) 

w 



This results explains why we have introduced the term coplanar. Note that aA ~ (3 
only when the corresponding invariant is much larger than mass 2 but the above 
expression is valid for all regimes and it is easily generalized to n emitters with the 
result that Q in a process 2 — > n any external charged leg % talks to all other charged 
legs, each time with a known scale and with a known total weight proportional to 



a (A\+...+A\)-l 



/r(a (A\ + ... + Ai)), 0< Xi <l (16) 



Note that each A has the appropriate sign, in/out, part/antp. Furthermore, is 
the number of pairs < ij > with i fixed. The IR exponent is given by 

2a fl+r 2 , 1 1 m? r , 
aA = — 1 2 111 -- 1 ' ITT = 71 Si 17 

7r [ 1 — r z r J |t| (1 — r) z 
For Bhabha scattering we will have the following combination: 

2 r i 

- A(s, m e ) - A(t, m e ) + A(u, m e ) = - In— ^--l, (18) 

71 



obtained as an exact result, not a guess. 



4.1 Conclusions for QED 



The structure-function language is still applicable but initial state structure func- 
tions evaluated for one scale is, quite obviously, not enough. In any process each 



1 A.Ballestrero, G.P. work in progress 
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external leg brings one structure function; since all charged legs talk to each other, 
each SF is not function of one ad hoc scale but all < ij > scales enter into SFj. The 
exact spectral-function is a convolution of SF 

£ P air<y>^ = [ ^K' ${K')E*f <ij> {K - K') , 

$(AT) = — — fd 4 x exp{iK-x + a (/-/«-,)} = 5{K) + 0(a). (19) 

(2 7TJ 4 J 

Furthermore, IR-finite reminders and virtual parts can be added according to the 
standard approach of reorganizing the perturbative expansion. 



5 Multi-loop calculations 

The third evolutionary development can be summarized as follows. A systematic 
approach to multi-loop calculations is needed. Although a lot can be done by means 
of approximations do we have a new, implementable algorithm for a full scale attack 
to the problem? 



There is a new idea around: the Tkachov theorem or generalized Bernstein 
theorem (see also Bardin's talk). 

V(x,d)l[Vr + \x) = B\{Vt{x). 

i i 

where V is a polynomial of vector of Feynman parameters) and 

di = d/df, B and all coefficients of V are polynomials of /i, and of the coefficients of 
Vi(x). 

The message is that we can do numerical calculations. However, one-loop |16| 
versus multi-loop is exactly as 2d Ising model versus 3d Ising model. For one-loop we 
have a universal master formula due to F. V. Tkachov [|15|] : 

V(x) = x t Hx + 2K t x + L, P = l- + d * 

2 + 1) 

X = K l H-\ B = L-K t H- 1 K. 
An example that everybody can do, the scalar 3-point function 

Co = t dx f dyV-^ix, y) = C d Q + C° + C° , 
Jo Jo 

V(x,y) = ax 2 + by 2 + cxy + dx + ey + f — i e 

C d o = 2 I 1 dx f X dyV(x,y)\nV(x,y)G 2 2 A- 2 , 
Jo Jo 
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G S = - J dx{V(x,x)\nV(x,x)[G 2 A- 2 (a x -a y ) + -A- 1 (a x -a y )G 11 ^i 1 
+ V(x, 0) In V(x, 0) [G 2 A- 2 a y + ^ A~ 1 a y G 12 A 2 1 

- 1/(1, x) In V(l, x) [G 2 A- 2 a x + ^ G^G^A' 1 + \g\A~ 2 + ^ A^a^A' 1 ] 
4 G^A-WA- + c + e) - 1 GaA-^xsAg 1 A + c + e) 

4 o o 

13 1 

+ — G 2 A~ 2 (3a + 6 + -c + 4d + 2e + 6/) + -A -1 a x ai x Ai x (a + 6 + c + d + e) 

J ^ 1 

1 1 3 3 
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4 o 4 4 
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1 1 3 3 

+ - A^a^A" x (a + d) - - A~ 1 a y G ll A- 1 (a + b + c+-d+-e) 

1 3 
+ -A- l a y G l2 A 2 \a + -d) 

+ - V(0, 0) In V (0, 0) [A-^^a^A- 1 - A^a^A" 1 + A^a^A" 1 

- J 1/(1, 0) In 1/(1, 0) [GaA-^gxAj 1 + A^a^A^ 1 + 4 A" 1 ^^ 1 
+ A-^^G^AJ 1 

+ J 1/(1, 1) In 1/(1, 1) [GaA-WAj 1 + G 2 A- 1 G 13 A3 1 - A^a^A^ 1 
+ A^a^as^A^ 1 

- A^a^GnA^ 1 + A^a x G 13 A 3 1 A~ 1 a y a lx A^ 1 + A^ayGnA^ 1 
Here, the Gram determinants are given by 

oja o J d + e d c + e 

a x = Idb - ce, a y = lea - etc, = — - — , a 2x = -, a 3x = - 

A = fG 2 - (bd 2 -cde + ae 2 ), 
(d + e) 2 



d 2 (r + e) 2 

— , A 2 = fG l2 --, A 3 = (a + d + f)G 13 -^L., 



Ai = /Gn- 

G 2 = 4a6 - c 2 , Gn = a + 6 + c, G 
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a, Gi 3 = 6, 



Another example is given by the IR-divergent Go- functions. One can prove that 
B = for 



a = b = f = m 2 , c = e = s — 2 m 2 , d = —2m 2 



As a consequence we have 
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d x y d y 
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This relation we use to write 
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Another example is the massless pentagon. Here (but also for the massive one) one 
can prove that 
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Now we raise the exponent and integrate by parts with B = 1/16 S12S23S34S45S51. 

In moving to multi-loop we recall that any diagram G with Nl legs and rii loops 
is representable as 
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r] is is the projection of line % along the loop s. T is a co-tree. Note that UV- singularities 
come from U so that, for finite diagrams, one should raise only the factor 

V = u\^m 2 iZi -Y. z i\ ~ Yl B n * ' Qi z ^ z i ■ 

i i ij 

There are several comments. For two-loops U is quadratic and V is cubic, so that we 
have to construct a Bernstein functional relation for a quintic or higher polynomial. 
Try to do it naively in FORM and you will be told 'Input expansion buffer overflow. 
Try something else.' 

For Nl = 3,4, we have V to a positive power and it is enough to regularize G 
with some Ks operation (S G G). For Nl > 5, we have V to a negative power, U 
to a positive power (no UV divergency), so we could raise a cubic. For self-energies 
at q 2 = it is easy. Consider, for example, rii = 2, Nl = 3, perform a projective 
transformation 

Zi = AiUi/^^AjUj, Ai — jj> 



mi 



to obtain V — — U quadratic and, 

V - (—^— — J—)- 

V X1J y „ o 12' 11/ "^22 

m\m\ mfm^ m^m^ mfm^ 

2 1 1 1 2 1 

2 2 22 22/~^~^ 22 

m\m2 m\m2 m\m^ m^m^ 

The crucial point being that we know how to raise a quadratic. Finally, Nl = 2 

is nice, a single polynomial of degree + 1. Of course, vacuum diagrams are also 

easy. At two-loop level we have 

G v 2 ac = (47t) £ ~ 4 r (e - 1) (K'H-'Ky 2 £ dx J* dy 



x 



£ — 2 i L £ - 1 



where C/ has been presented above. The same is, obviously true, for tadpoles. What 
we are still missing is a real multi-loop example, 

Additional problems are represented by the fact that V is an highly incomplete 
polynomial, few non-zero coefficients. V can be found via a direct study of linear 
systems of n equations (n ^> 100) in m variables and very often the matrix of co- 
efficients has rank < n, so the system is generally impossible and, we need to go to 
higher n, m, i.e. beyond n min = min{n < m} (try one-loop master formulae with 
K = L = 0!) and your algebraic manipulation will exhaust your computer before you 
find a solution. 
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A simple counting shows the following situation: given 
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four variables and a cubic V require n eq = 126(330, 715) and n Yai = 155(415, 871) for 
n = 2(1,0) and first(second, third) derivatives. Simple examples already show that, 
for realistic polynomials V, one has to go beyond second order in derivatives. 

The relation between B = and Landau singularities is still to be examined, 
although B = contains all singularities of G, physical and spurious ones. L — 
K l B^ x K (or B in general) is not a Gram determinant but its spurious zeros, if any, 
are a vexation, as much as those of Gram determinants in the standard reduction 
procedure. For Co (— m 2 , —ml, —s; mb, M, m a ) we obtain 
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